Emergence of Chiral Magnetism in Spinor Bose-Einstein Condensate 

with Rashba Coupling 
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Hydrodynamic theory of the spinor BEC condensate with Rashba spin-orbit coupling is presented. 
A close mathematical analogy of the Rashba-BEC model to the recently developed theory of chiral 
magnetism is found. Hydrodynamic equations for mass density, superfluid velocity, and the local 
magnetization are derived. The mass current is shown to contain an extra term proportional to the 
magnetization direction, as a result of the Rashba coupling. Elementary excitations around the two 
known ground states of the Rashba-BEC Hamiltonian, the plane-wave and the stripe states, are 
worked out in the hydrodynamic framework, highlighting the cross-coupling of spin and superfiow 
velocity excitations due to the Rashba term. 
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A recent excitement in the cold atom physics is the 
newly found ability to engineer spin-orbit-type interac- 
tions among the spinor 87 Rb atoms [1]. The effort is 
part of a broader theme to create synthetic gauge field 
environment for cold atoms [2] and currently forms one 
of the most exciting branches of cold atom research. Sev- 
eral theoretical papers appeared dealing with the possible 
phases of spinor Bose-Einstein condensate (BEC) in the 
presence of spin-orbit interaction of the type ~ vp T F ■ pvf», 
for an appropriate spin operator F, the spinor ^, and mo- 
mentum p = —iTiW [3, 4]. Due to the similarity of this 
term to the Rashba coupling ~ cr x Py — &yPx (which can 
be mapped to a x p x + <J y p y with suitable spin rotation) in 
solid state physics [5], spinor BECs containing such in- 
teraction will be called Rashba-BECs. The Rashba cou- 
pling contains one space gradient while the usual kinetic 
energy has two. The competition between the two ener- 
gies results in a length scale, as manifested for instance in 
the striped superfluid state in certain interaction param- 
eter regimes of the Rashba-BEC Hamiltonian [3]. The 
influence of confining trap on the Rashba-BEC was also 
considered [6]. In another vein, the rotation effect in a 
harmonically trapped Rashba-BEC was studied in sev- 
eral papers [7]. Fast-rotating limit of the Rashba-BEC 
is expected to yield non-Abelian Landau level structures 
with the exotic possibility of a many-body state charac- 
terized by fractional statistics [8]. 

Most theoretical works so far on Rashba-BEC have 
been numerical studies of the ground states. Recently, 
an interesting possibility of fractionalized vortex excita- 
tion in the striped superfluid state and associated ther- 
modynamics of Kosterlitz-Thouless phase transition were 
considered in Ref. [9]. Still, the general characterization 
of the Rashba-BEC system in hydrodynamic variables 
and identification of proper conservation laws are lack- 
ing. In Ref. [10] the authors studied fluctuations of the 
so-called plane-wave state within the framework of the 
time-dependent Gross-Pitaevskii (CP) equation directly, 
while the discussion of more general states in spin-orbit 



coupled BECs could be found in Ref. [11]. In this paper, 
in the tradition of formulating superfluid phenomena in 
hydrodynamic variables [12], we derive conservation laws 
of mass density (p) and mass current (J) for Rashba- 
BEC system. In addition, Euler equation for the super- 
flow velocity (v) and Landau-Lifshitz equation for the 
spin vector (n) dynamics are obtained. These are gen- 
eralizations of the hydrodynamic formulation carried out 
for non- Rashba spinor BEC case in Refs. [13, 14]. As will 
be shown, spin-orbit coupling introduces several new fea- 
tures missed in previous works such as the rigorous math- 
ematical mapping to chiral magnets and the modification 
in the meaning of the current density. The set of equa- 
tions are then applied to study coupled fluctuation dy- 
namics of the hydrodynamic coordinates for plane-wave 
and striped superfluid states of Rashba-BEC Hamilto- 
nian. Hydrodynamic theory provides a useful alternative 
to the Bogoliubov approach to study fluctuations [15] and 
our calculations illustrate the coupled nature of spin and 
superfiow velocity in the excitation spectrum. 

We begin by writing down the single particle Hamilto- 
nian of Rashba-BEC system as 



Ho = l /[(p + «F)*]t-[(P + «F)*]» 



(1) 



where k describes the strength of Rashba spin-orbit cou- 
pling. We will assume h = m = 1 throughout the paper. 
The task at hand is to transform the Rashba-BEC Hamil- 
tonian Hq using hydrodynamic variables. For spin-1/2 
case in particular, we have F = <r (2 x 2 Pauli matrices) 
and one can decompose the two-component spinor field 
* as [13, 14] 



i/jz, 



cos(6»/2)e-^/ 2 
sin(6>/2)e^/ 2 



(2) 



where ip = ■s/pe^ with £ being the overall phase of 
the spinor condensate. The CP 1 field z appears in the 
theory of nonlinear a- model (NLerM). The relevant hy- 
drodynamic variables are defined in terms of the spinor 
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field as p = *tvj/ and n = z+crz = ^cr^/p = 
(sin 9 cos (j>, sin 9 sin <fi, cos 0) . Superfiow velocity is de- 
fined asv = a c -a s , where the two vector potentials are 
a c = — ie -t ^Ve*£ = V£ and a s = iz^Vz. Certain identi- 
ties can be derived relating terms in the GP Hamiltonian 
to the hydrodynamic expressions as [16, 17]: 

(ptt)t • (p*) = p(± * n • 9 ** + v ') + ( V VP) 2 , 

i 

l^er • (p*) + h.c. = p(n ■(Vxn)|2n-vj. (3) 

The Rashba-BEC Hamiltonian (1) is accordingly trans- 
formed to H = Hi + H 2 + H 3 where 

Hl = f [\ 9in ' d i n + Kn ' (V x n) 



^ 2 = i(vvp) 2 , 
f 3 = \ (v+H 2 - 



(4) 



A similar exercise for spin-1 operator F and ferromag- 
netic spinor wave function [18] yields the same hydrody- 
namic Hamiltonian as Eq. (4), only with 1/4 in Hi being 
replaced by 1/2. The analysis carried out below therefore 
applies equally well to both two- and three-component 
spinor BECs with Rashba coupling. 

The key new element in the spin dynamics (Hi) in- 
troduced by nonzero k is n • (V x n). In solid state 
magnetism, such a term is known as the Dzyaloshinskii- 
Moriya (DM) interaction. Taken together with the 
NL<jM term, the ground state of the spin Hamiltonian 
(1 /4) J2i di n ■ din + nn ■ ( V x n) is the well-known spiral 
spin phase with the modulation wave vector |q| = 2k, and 
the spins rotating in the plane perpendicular to q. Here 
we have 1/4 for spin stiffness J in magnetism, and k as 
the ratio of DM energy D over J. As we will show shortly, 
the stripe state [3] carries precisely such spin structure. 
In the incompressible limit, assuming v + nn = to min- 
imize H$, we find the perfect match of the Rashba-BEC 
Hamiltonian to that used in the study of spiral mag- 
netism [17]. In the NLctM+DM spin model, a textured 
spin phase called the Skyrmion crystal was identified in 
the presence of weak magnetic field or spin anisotropy 
(See Ref. [17] and articles cited therein). We believe it 
is no coincidence that a very similar Skyrmion lattice 
was found in two recent theoretical studies of the phase 
diagrams of Rashba-BEC Hamiltonian in a trap [6] . 

Next we present the dynamical equations of motion 
for p, v, and n, based on the Hamiltonian H — Ho + Hi 
where Hq is given in Eq. (1) and Hj, the interaction, is 
given by 



Hi 



[9c{Pt + Pif ~ 9s{Pt ~ Pi) 2 ] 



(5) 



for spin-1/2 case. Here p a (a =t,|) refers to the density 
of the spin-CT component. Note that g c = (g^ + g^i)N/2, 



g s — (g-fi — g^)N /2, with and g^i being the intra- 
and inter-component s-wave interactions, respectively, 
and N is the total particle number. It is assumed that 
0ft = 9U an d 5fi — 5j.f ■ We state the equations first, 
and discuss briefly methods of derivation later. 



(1) Mass conservation equation: 

dtp + V ■ [p(v + sn)] = 0. 



(6) 



This is the generalization of the usual mass conservation 
equation to the case with nonzero Rashba coupling 
parameter k. The combined expression p(v + nn) = J 
serves as the mass current density. The relationship can 
be expected from taking a functional derivative of the 
hydrodynamic Hamiltonian Ho given in Eq. (4) with 
respect to v: 5Hq/5v = SH^/Sv = p(v + nn). Such 
definition of the current density, albeit odd at first sight, 
was anticipated on symmetry grounds by Ambegaokar et 
al., who constructed the Ginzburg-Landau functional for 
the A- phase of superfluid 3 He [19]. They argued that the 
definition of the current density in general should consist 
of three types of terms, the first two of which would 
reduce to pv and npn in the case of isotropic superfluid 
density ptj — pSij. The third term, proportional to 
the magnetization current V x n, does not exist in our 
theory but might appear when a more general Hamilto- 
nian is employed. It was mentioned [19] that n (axial 
vector) can appear in the definition of J (polar vector) 
only if the medium breaks the inversion symmetry. The 
Rashba term, being odd under V — > — V, breaks such 
symmetry. Additionally, the newly defined J satisfies 
the requirement of being gauge- invariant, if we view the 
Rashba coupling as a non-Abelian gauge field. 



(2) Euler equation: 

D t w = v x b + e Vp, 

where 

p= -^(9 4 n) 2 + Kn- v+-n- (Vxn) 



(7) 



V 2 ^ 



+ 9cP ~ 9spn z 



2VP 

The material derivative D t = d t + v • V appears in 
the above. Compared with the case in the absence of 
Rashba coupling [13, 14], only the explicit expression 
of the quantum pressure p is modified by k. The 
internal electric (e) and magnetic (b) fields arise from 
spatial and temporal fluctuations of the magnetiza- 
tion in the usual way [14]: a — — • (9,n x ftn), 
h = -\e ijk Yi ■ (djn x d k n). 

(3) Landau- Lifshitz equation: 
pD t n = nx ( l-d i (pd l n)-2Kpv- kV x (pn) + 2g s p 2 n z 



(8) 
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where n z = (0,0, n z ). Compared to the earlier expres- 
sion [14], two new terms appear on the r.h.s. due to the 
nonzero re as in Euler equation. 

The three equations (6), (7), and (8) together describe 
the dynamics of Rashba-BEC system in hydrodynamic 
variables. For their derivations one can substitute Eq. (2) 
into the time-dependent GP equation 



^(-iV+K*) 2 +g c (&*)-g s (&<T z *)a 2 



and follow the projection strategies outlined in Ref. [14]. 
Alternatively, one can start from the hydrodynamic La- 
grangian 

C = p[-dt£ + iz%x] - H[p, v, n] + A(zt z - 1), 

the constraint z^z = 1 being implemented by the La- 
grange multiplier field A. The variation of the La- 
grangian C with respect to £, p, and zt, respec- 
tively, yields the mass continuity, Euler equation, 
and spin dynamics as shown above. The Lagrange 
multiplier A term is eliminated in the final equa- 
tion through the projection with ztr, where ztr = 
(cos(6»TR/2)e-^ TR / 2 ,sin(6»TR/2)e^™/ 2 ) T with 9 TR = 
7r — 6 and ^>tr = <f> — 7r, is the time reversed version 
of z satisfying z^ R z = 0. 

As discussed earlier, the mass conservation law shown 
in Eq. (6) suggests the definition of the mass current 
density J = p(v + /en). In a closed system, the current 
density J is expected to obey its own continuity equation 
d t Ji + djUij = with an appropriate momentum flux 
tensor Lly . After some lengthy algebra, we arrive at such 
expression as: 



17 — P 



(9) 
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Pij = ~ (V 2 p) Sij + pviVj + ^din ■ dj-a 

Qij = 2np5ij + p{viUj + mvj) 

-^([nx d,n} 3 + [nx d 3 n]^. 

We have separated the re-dependent part as in the 
above, while the re- independent part is denoted as P^. 
The last term in Qij is the well-known spin current in 
spiral magnetism [20]. It ought to be emphasized that 
the continuity equation applies for J = p(v + ren), but 
not for pv alone. Such observation strengthens our earlier 
claim about the proper definition of mass current in the 
Rashba-BEC medium. 

Hydrodynamic equations such as derived above can be 
applied to study small fluctuations of the known ground 
states. The three coupled equations (6), (7) and (8) have 
to be solved simultaneously, while also obeying the con- 
straints implied by the Mermin-Ho relation [21], 



Most small fluctuation theories such as carried out in 
Ref. [15] deal with the time-dependent GP equation in 
the Bogoliubov approach without explicitly introducing 
the Mermin-Ho constraint. Here we establish that the 
hydrodynamic analysis, carried out in a way to consis- 
tently take care of the Mermin-Ho constraint, yield iden- 
tical results as the standard Bogoliubov approach. The 
great advantage of hydrodynamic analysis, on the other 
hand, is that it is much easier to identify the nature of the 
given excitation mode as either density, spin, or super- 
fluid velocity oscillations, or some mixture thereof. Be- 
low we examine two known ground states of the Rashba- 
BEC Hamiltonian H = Hq + Hj. The first one, called 
the plane- wave state, is obtained for g s < and reads 
*o = (l/\/2)e tK:E (l,-l) T [3]. The other stripe state ex- 
ists for g s > 0, and reads — (cos rex, —i sin rex) T [3]. 
The two ground states become degenerate at the critical 
point g s = 0, with the energy E = — n 2 /2 + g c . In typical 
experimental situations we always have g c ^> \g s \, hence 
we take g s = in the following to simplify the analysis. 

The plane- wave state is characterized by po = 1, 
vq = (re, 0,0), and no = (—1,0,0). It is a state with 
ferromagnetic spin order no and nonzero superflow vo, 
which however combine to produce zero mass current 
Jo = vo + reno = 0. Small fluctuations around each 
variable is parameterized as p = po + X, v = vq + v, 
n = n +T] y e y +r] z e z , and inserted into the hydrodynamic 
equations (6), (7) and (8) to first order in Xt^tVviVz- 
Assuming the two-dimensional fluid, fluctuation in the 
z-dircction will be ignored. The effective electric and 
magnetic fields do not contribute to the Euler equation 
in the linear analysis. Furthermore, Mermin-Ho relation 
as applied to the plane-wave state yields the constraint 
V x v = 0, which is solved by writing v = V/. In k- 
space, the scalar functions % an d /> an d two components 
of the spin fluctuation rj y and rj z obey a four-dimensional 
matrix equation of motion, M(x, f, Vy> Vz) T — 0, where 
M. equals 



V 



-uj i(k%+ky) nk y 
--ig c -co 





2 2 

2nk x — ijj —i^- — 2K 2 i 



J 



i 



(V x v)j = -j-e^fen • {dj-a x d k n). 



(10) 



Diagonalizing the matrix, one obtains the dependence 
of excitation mode uj on k as shown in Figs. 1(a) and 
1(b). For the parameters we choose, we can observe the 
branches of one shifted mode in k x direction and one 
gapped mode in k y direction. These features coincide 
with the ones obtained in Ref. [10], which was based 
on the analysis of time-dependent GP equation in the 
Bogoliubov approach. Furthermore, we label in the 
dispersions the two components which contribute the 
most and the second most to each eigenvector. For 
instance, the shifted modes in Fig. 1(a) come from the 
spin fluctuation, while the linear sound modes are due 



4 



to fluctuations of the total density and the overall phase. 
Similar analysis applies to Fig. 1(b). Hydrodynamic 
theory indeed enables us to directly identify different 
modes more easily than Bogoliubov approach. 
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FIG. 1. (color online) Excitation modes of plane-wave (up- 
per panels) and stripe (lower panels) state along the k x (left 
column, k y = 0) and k y (right column, k x = 0) directions. 
Rashba strength k = 1 sets the unit for the z-axis and g c — 10. 
Only the four gapless branches are shown in the lower panels 
for the stripe state. We choose g s — for all cases. Shown 
inside the parentheses are the two components contributing 
the most and the second most to each eigenvector. 



pling of the k-modes displaced by ±q is a consequence of 
the ground state being a modulated state, thus no way 
to write the excitation equations into the 4x4 matrix 
form as in the plane wave case. Numerical diagonaliza- 
tion yields the spectrum shown in Fig. 1(c) and 1(d). 
Analogous "helimagnon band" was identified in a spiral 
magnet MnSi [22]. For the parameters we choose, four 
linear sound modes with mode contributions as labeled 
in Fig. 1(c) and 1(d) exist near k x = k y = 0, which is 
different from the plane-wave case. Bogoliubov analysis 
of the GP equation yields identical spectra for the stripe 
state as well. 

In summary, we have derived the mathematical map- 
ping of the Rashba spin-orbit-coupled BEC Hamiltonian 
in terms of hydrodynamic variables. The mapping high- 
lights the connection of the present system to the chiral 
magnets. The hydrodynamic equation for mass conserva- 
tion, as well as Euler and Landau-Lifshitz equations for 
superfluid velocity and spin dynamics, respectively, are 
derived. The mass current contains an extra magnetiza- 
tion term due to the Rashba coupling. As applications 
of hydrodynamic approach, we studied the elementary 
excitation modes of the plane-wave and stripe states of 
Rashba-BEC, featuring the different behaviors of fluctu- 
ations around the two degenerate ground states and the 
coupled behavior of spin and superflow velocity in the 
presence of Rashba coupling. 

H. J. H. is supported by NRF grant (No. 2011- 
0015631). Useful discussion with Hui Zhai at the early 
stage of research is gratefully acknowledged. 



The stripe solution yields po = 1, vo = 0, and 
n o = (0, — sm9, cos 9), where 9 = 2kx. The spin order 
is a spiral with the plane of spin rotation orthogonal to 
the propagation vector q = 2ne x . Only the fluctuations 
of n need to be parameterized differently from the plane- 
wave case, as n = n + J7iei + rj2e 2 , where ei = (1, 0, 0), 
e2 = (0, cos 9, sin 8), and no form the three orthogonal 
axes in the rotating Frenet-Serret frame [13]. Contrary to 
the plane-wave case, we have a nonzero effective electric 
field contribution to the Euler equation as e x — —nrji and 
correspondingly, the Mermin-Ho constraint in the case of 
stripe state is solved by v x -\-Kt)\ = d x f, v y = d y f. Insert- 
ing them back into the hydrodynamic equations yields 



-ujx + 


1 u - 
2 Kk vX 


- + ik 2 f - 


K + 

h Q kyVi 


= 0, 


- ig c )x 


-uf + 


2 K ^yf 


«. + 


= 0, 




- WTIi - 


i , - 

-KKyT] 1 


i 1 2 
+ - 2 ^m 


= 0, 


- 2iK 2 r]i 


- UJTI2 " 


h Kk x rj 2 - 


i , - 


= 0, 



where — m(k — q) ± m(k + q) (m = Xi /? ViiVz) an d 
un-superscripted variables are at momentum k. Cou- 
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